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Abstract 



We show the theoretical and experimental combination of acoustic and 
optical methods for the in situ quantitative evaluation of the density, the vis- 
cosity and the thickness of soft layers adsorbed on chemically tailored metal 
surfaces. For the highest sensitivity and an operation in liquids, a Love mode 
surface acoustic wave (SAW) sensor with a hydrophobized gold coated sens- 
ing area is the acoustic method, while surface plasmon resonance (SPR) on 
the same gold surface as the optical method is monitored simultaneously in a 
single set-up for the real-time and label-free measurement of the parameters 
of adsorbed soft layers, which means for layers with a predominant viscous 
behavior. A general mathematical modeling in equivalent viscoelastic trans- 
mission lines is presented to determine the correlation between experimental 
SAW signal shifts and the waveguide structure including the presence of 
the adsorbed layer and the supporting liquid from which it segregates. A 
methodology is presented to identify from SAW and SPR simulations the 
parameters representatives of the soft layer. During the absorption of a soft 
layer, thickness or viscosity changes are observed in the experimental ratio of 
the SAW signal attenuation to the SAW signal phase and are correlated with 
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the theoretical model. As application example, the simulation method is 
applied to study the thermal behavior of physisorbed PNIPAAm, a polymer 
whose conformation is sensitive to temperature, under a cycling variation 
of temperature between 20 and 40 °C. Under the assumption of the bulk 
density and the bulk refractive index of PNIPAAm, thickness and viscosity 
of the film are obtained from simulations; the viscosity is correlated to the 
solvent content of the physisorbed layer. 
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1 Introduction 

Sorption processes at the solid/liquid interface by which (bio)molecules bound 
to material surfaces are of interest for biosensors, biomaterials, materials 
characterization and surface science. Understanding the three-dimensional 
organization (including density, water content and thickness) of the resulting 
sorbed film and its evolution during the process is crucial for many applica- 
tions in these domains. For biosensors, more specifically, there is a need to 
monitor the response in real-time and to be able to distinguish contributions 
coming from the dry sorbed mass, which is the physical criterion for esti- 
mating sensitivity, and those that should be attributed to effects intimately 
associated to the film organization, like sorbent-bound water and viscosity 
for example. While a wide variety of methods can qualitatively detect the 
formation of sorbed films, almost none of them, alone, is quantitative and 
able to reveal the film organization. Furthermore, only a few techniques 
can monitor the sorption process in real-time. Scanning probe microscopies 
might fulfill all these requirements, mainly for submonolayers sorbed on ul- 
tra smooth surfaces [T]. Neutron reflectivity XPS j^j, mass spectroscopy 
jl] and radiolabeling j3] are quantitative techniques able to directly measure 
the dry sorbed amount, nevertheless compatible to very different film orga- 
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nizations. Of all the direct detection (i.e. label-less) techniques, we have 
identified acoustic and optical methods as being the only ones fulfilling two 
fundamental criteria of our measurements: time resolved and in situ (liquid 
phase) measurement of the physical properties of the adsorbed layer. 

Various methods of direct detection of biochemical layers have been devel- 
oped, either based on the disturbance of an acoustic wave [Bj (quartz crystal 
microbalance [7] - QCM - and surface acoustic wave devices jHl HH HH 
H2J EH Ulj - SAW) or of an evanescent electromagnetic wave (optical waveg- 
uide sensors [13 EH and surface plasmon resonance [T71 [TH] - SPR). While 
each one of these transducers provides reliable qualitative curves during the 
adsorption of a soft layer on their functionalized surface, such as proteins 
for example, extraction of quantitative physical parameters such as optical 
index, density, viscosity or water content requires modeling of the adsorbed 
layers [T^]. The modeling includes multiple parameters which must be iden- 
tified simultaneously: hence the need for the combination of acoustic and 
optical detection methods in a single instrument (221 ED E3 120] • 

We here use a combination of Love mode surface acoustic wave (SAW) 



5 



device and surface plasmon resonance (SPR) and present a methodology to 
extract out of experimental data and simulations a quantitative value for 
the density, the viscosity and the thickness of adsorbed soft layers from a 
liquid phase. The viscosity and the density will be further correlated to the 
water content in the adsorbed layer. As a showcase of the performances 
of our methodology, the detection of the physisorption of poly(N-isopropyl 
acrylamide) - PNIPAAm - is presented. PNIPAAm is a thermally sensitive 
polymer with a conformal transition investigated previously by SPR and 
QCM techniques !2Uj. The chemical structure of PNIPAAm is reported 
in Fig. [TJ Our approach is presented as a sequel to the study by Plunkett 
et al. [201 where the properties of PNIPAAm physisorbtion to a hydropho- 
bic surface was monitored by a combination of SPR and QCM, and the 
subsequent temperature-dependent conformation transition was monitored 
by QCM alone by measuring both the resonance frequency and the qual- 
ity factor of the resonator. In the present study, Love mode SAW devices 
were chosen for a high mass sensitivity and a compatibility with measure- 
ments in liquid media [23 I2E1 ESI I2H]- Being based on the propagation of a 
shear-horizontal acoustic wave, their interaction with the surrounding liquid 
is minimal although the bulk properties of the solvent influence the acoustic 
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wave propagation and must be considered in the analysis of the sensing device 
[23 HH] and a model for the analysis of the structure, including the adsorbed 
layer and the liquid, is presented in the theoretical section (section |2J). While 
both SAW and SPR sensors display theoretically similar sensitivities and de- 
tection limits, acoustic sensors are typically sensitive to temperature fluctua- 
tions unless optimized to achieve a zero temperature coefficient of frequency 
(TCF). In a normal operation, care must be taken to properly control the 
environment around the instrument to keep a stable temperature. However, 
in the present case, the temperature will be voluntary modified and moni- 
tored to modify the characteristics of the adsorbed PNIPAAm. This specific 
setting requires a preliminary calibration curve of the acoustic response as 
a function of the temperature. The analysis of the layer will assume that a 
certain amount of the solvent - water in the present case - is contained in 
the layer. A schematic representation of the experimental combined set-up 
used for this study is given in Fig. El 



2 Theory 

Love mode surface acoustic waveguide. Acoustic devices are sensitive 
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to mechanical perturbations, temperature changes and electrical perturba- 
tions such as pH changes, conductivity and dielectric permittivity of added 
materials. If the sensing area is electrically shielded by a grounded con- 
ducting layer, these electrical effects can be neglected, letting the device be 
sensitive to mainly mechanical perturbation. For the model, we will assume 
that a composite material segregates from the bulk solution on the waveg- 
uide surface as an hypothetical perfectly flat layer of given characteristics to 
be determined. Acoustic waves with a shear-horizontal polarization are best 
suited for sensing applications in liquids: the mechanical displacement at the 
surface probes within a certain volume the density, the shear stiffness and the 
dynamic viscosity of the composite material that has a certain thickness and 
is adsorbed on the surface, thus resulting in a mechanical interaction between 
the layer and the acoustic device. The Love mode acoustic structure sup- 
ports the guided propagation of a shear-horizontal wave in a stack of guiding 
layers deposited atop a piezoelectric crystal substrate. For a large velocity 
difference between the guiding layers and the substrate, the acoustic energy 
is trapped closely to the sensing area and results in the highest sensitivity 
of all similar acoustic structures. The acoustic velocity and the attenuation 
of the acoustic wave in the stacked structure is electrically probed thanks to 
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interdigitated electrodes patterned on the piezoelectric material. The phase 
at a given frequency within the bandpass range of the delay line and the 
insertion loss of the acoustic device measured electrically, for instance with a 
network analyzer, are representative of the actual propagation characteristics 
of the acoustic wave along the delay line and, in particular, along the sensing 
path. 

Equivalent transmission line model. It is known that the propagation 
of acoustic waves in guiding structures can be obtained by equivalent trans- 
mission lines pIT) 132] • An implementation of this equivalent model is given 
here to analyze the propagation characteristics of the guided wave in the 
Love mode structure, including the presence of the surface adsorbed com- 
posite layer and the liquid from which this layer segregates. At first, each 
material is equivalently represented by a transmission line with a wavevector 
k and a characteristic impedance Z such that the shear stress T(r) and shear 
particle velocity v(r) fields at a given point f in the material are given by 
the superposition of incident and reflected acoustic waves: 



T(f) 
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where T + and T_ are arbitrary values for the intensity of the incident and 
reflected waves. The equivalent representation of a viscoelastic layer in this 
model is given in Fig. EJ which consists in a series inductance L shunted by 
a capacitance C placed in series with a conductance G. The values of the 
equivalent parameters are given in mechanical terms of the layer, i.e. its 
density p, its shear stiffness p and its dynamic viscosity r\ by the following 
correspondence [32] 

L = p 

C = - 

P 

G = - 
V 

from which we derive an analytical expression of the wavenumber k 
and of the characteristic impedance Z: 

k - - 
k ~ V 

p + iurj 
V 

where u = 2nf is the angular frequency, while V is the complex acoustic 
velocity defined by 

lp+i,,'i) 
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For waveguides, the wavenumber k is decomposed in a resonance wavenumber 
k z and a propagation wavenumber k x such that k 2 = k^ + k 2 . The projection 
into both directions, propagation and resonance, is expressed in each layer 
by a complex coupling angle if such that 

kj X = kjsimpj (1) 

kj z = kj cos (pj (2) 

Zj X = Zj sin (pj (3) 

Zj Z = Zj cos tfj (4) 

where the subscript j refers to a constituent layer. The transverse resonance 
principle^n\ is applied to determine the propagation wavenumber: according 
to this principle, the coupling angles are related between each others such 
that the propagation wavenumbers kj X are all equal. For a simpler notation, 
the kj X being all identical are noted k x . 



To determine the propagation characteristics in the waveguide structure, 
the last step is to determine the resonance wavenumber. To that purpose, an 
equivalent circuit of the physical structure represented in Fig. 0] is built by 
connecting in series each equivalent transmission lines of the layers stacked in 
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the direction of resonance, as depicted in Fig. |5J The series connection auto- 
matically solves the boundary conditions associated to the continuity of the 
mechanical stress and strain at materials interfaces. As a consequence, semi- 
infinite layers (the substrate S and the fluid F) are terminated by a matching 
load that accounts for the endless side while the other layers (the guiding lay- 
ers Gj and the composite layers Cj) have finite dimensions, respectively tj 
and hj. One must notice that for a clearer representation, only one guiding 
layer and one composite layer have been schematically depicted although the 
model is easily extended to a larger amount of layers. Finally, the resonance 
is found by measuring the impedance at a plane T located at an arbitrary 
position in the structure, for instance at the substrate-coating interface as 
sketched in Fig. the resonance is found for values of the complex coupling 
angle that result in opposite impedances seen in both directions at the ref- 
erence plane T. These impedances are determined by applying a standard 
formula from the transmission line theory to calculate the input impedance 
of a transmission lines of characteristic impedance Z c h ar and length L loaded 
at its end with an impedance Zi oad : 

2 _ g Zioad + Z char tanh (k z L) 
Zchar + Zi oad tanh (k z L) 

A solution of the resonance wavenumber is translated to the propagation 
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wavenumber thanks to Eqs. (|TJ and (J2J). The real and imaginary parts of a 
determined solution k x are related to the attenuation and the phase velocity 
of a guided mode, respectively. 

Electrical signal. For the analysis of the acoustic signal, variations in 
amplitude and in phase are monitored during modifications of the mechanical 
characteristics of the composite layer or in its thickness. Because of such 
variations, the wavenumber solution changes to k' x . The induced variations 
of the acoustic signal are measured through the transducers and interpreted 
as an electric signal attenuation per unit of length A A, in dB/m, given by 



and to an electrical delay phase shift per unit of length A0, in 1/m, given by 

A(j) = lm(k x ) -lm(k' x ). 

The quantities to compare with experimental values are obtained by multi- 
plying the theoretical shifts by the length of the sensing path D. 



AA = 20 log 




Rigid and viscous interactions. The ratio of the amplitude shift to the 
phase shift is a significant indicator of rigid or viscous interactions between 
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the layer and the vibrating surface and will be used to discriminate the re- 
laxation period of the adlayer. For any added layer on the sensing surface, 
four parameters must be evaluated from the measurement (density, shear 
stiffness, viscosity and thickness) while the SAW device gives only two in- 
dependent measurements linked to the real and the imaginary part of the 
acoustic wavevector measured electrically by the transducers (insertion loss 
and phase); therefore, an assumption must be made to solve some complexity 
of the analysis to evaluate the unknown parameters. To that purpose, we 
assume that the composite layer and the surrounding liquid are Newtonian 
fluids, meaning that the shear stiffness in both layers, \ic and fip, is discarded 
in regard to the viscous components ujtjc and ur] F . This working hypothesis 
has been previously described in the literature [33J through the relaxation 
time t = rj/fi, and a comparison of the relaxation time to the excitation 
frequency uj driving the mechanical motion of the waveguide surface. The 
product of the relaxation time by the frequency indicates if the adsorbed 
film moves in phase (rigid) or out of phase (viscous) with the acoustically 
driven surface. When the relaxation time is high, i.e. when urj 3> //, the 
film is viscous, meaning that acoustic loss is a predominant mechanism of 
energy dissipation in the film, which is also a mechanism driving an impor- 
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tant attenuation of the acoustic wave as it propagates along the perturbed 
part. At the opposite, when the relaxation time is low and lot] <C /i, the film 
is considered as rigid, involving an entrapment of the acoustic energy in the 
film which participates effectively in the guiding of the acoustic wave. 

For viscous films, it has been demonstrated |HH that the product of the 
density by the viscosity and the film thickness play a role in the modification 
of the acoustic velocity and the attenuation of the signal; for rigid films, only 
the density and the film thickness play a role. A discrimination between 
both effects can be obtained through the ratio of the SAW signal AA/A<p. 
The mechanical behavior of the film with respect to the excitation frequency 
is determined by the part, real or imaginary, of the propagation wavevector 
k x that is predominant. For the operating frequencies of Love mode devices, 
typically comprised between 80 MHz and 430 MHz, a film is seen as a viscous 
layer till shear stiffness up to 10 5 Pa. Above this threshold, the film supports 
the propagation of the guided Love mode rather than its dissipation, result- 
ing in a low ratio AA/A<p and leading to the influence of the density of the 
film only rather than the product of the density by the viscosity. The value 
of the ratio is obtained from simulations for each specific Love mode device. 
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Penetration depth and film thickness. The penetration depth of shear 
waves in fluids 5 is an indication of the thickness probed by the acoustic 
method. The exact value of 5 for surface acoustic waveguides is determined 
by the real part of the wavevector in the viscous layer as determined by the 
transmission line equivalent model. Its value is close to the one obtained 
from the well-known relation [33] 



from which it slightly differs due to the coupling angle (pc in the viscous 
layer. Because of a usual large velocity difference between the rigid guiding 
layers and the viscous layers, the real part of the coupling angles (pc and 
ifp are almost equal to zero and, therefore, the penetration depth can be 
identified with the expression given by Eq. The acoustic wave senses 
the viscous layer with an exponential decay extending to a thickness about 
3(5, after which the acoustic field decays below 10% of its initial value and is 
neglected. This effect has an impact on the ratio AA/A<j) as plotted in Fig. EI 
where the ratio is plotted versus the thickness of the viscous film normalized 
to the penetration depth as simulated for different possible values of r/c- The 
ratio stabilises for film thicknesses above ~ 35, which is about 150 nm at a 




(5) 
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typical acoustic frequency of 123.5 MHz for a layer with rjc = 1 cP (same 
as water). When the layer is thin (< 0.15), the ratio becomes similar to the 
one obtained with a rigid layer; consequently, a small experimental ratio is 
not sufficient to conclude that the probed layer behaves as a viscous or a 
rigid layer. On the other hand, a high ratio is an unquestionable indicator 
of a viscous interaction, where the hypothesis of a soft layer can be assumed 
(corjc 3> A*c)- The estimation of the shear stiffness of the adsorbed film is 
removed from the film analysis while the density, the viscosity and the thick- 
ness remain to be determined. These three values can not be estimated from 
the two information (insertion loss and phase) resulting from a single SAW 
measurement. Consequently, an additional information about any of these 
three values is requested with a measurement technique compatible with sur- 
face acoustic waves. A first possibility is a thickness estimate by AFM 
A second possibility is an optical measurement of the film such as ellipsom- 
etry or surface plasmon resonance. This latter possibility has been chosen 
and is described here after. 

Combined technique SAW/SPR. We have previously demonstrated the 
integration of a SAW device with surface plasmon resonance (SPR) in a single 
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instrument for the simultaneous and real-time measurement of rigid adsorbed 
layers jHHl- The SPR signal is a measurement of the coupling plasmon angle as 
a function of variations in refractive index n and thickness h of the adsorbed 
layer. A single SPR measurement is therefore not sufficient to determine 
independently these two parameters. Our previous analysis applied to rigid 
layers, therefore the SAW shift was assumed to be proportional to the density 
and thickness of the layer with respect to a mass sensitivity calibration of 
the SAW device by copper electrodeposition. In the rigid layer assumption, 
the acoustic phase difference combined to the SPR signal is sufficient to 
estimate independently the density and the thickness of the layer under the 
assumption of a linear relationship between the variation of the density of 
the layer and its refractive index variation. This assumption, which we will 
assume that still hold in the present case, considers that both the density 
and the refractive index of the composite layer are directly correlated and 
vary between the bulk values of the component segregating from the fluid on 
the sensing surface alone and the bulk values of the fluid itself as a function 
of the quantity of fluid contained in the composite layer. If we denote by the 
term x the relative amount of fluid in the composite layer, meaning that x 
varies between and 1, then the density pc of the composite layer is given 
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by 



pc = (1 - x)p L + xp F 



(6) 



and its refractive index nc by 



nc — (1 — ^)riL + xra^ 



(7) 



where the subscript F refers to the semi-infinite fluid layer, C to the compos- 
ite surface adsorbed layer, and L to the bulk characteristics of the compo- 
nent segregating from the fluid. By eliminating the unknown term x between 
Eqs. El and a direct relation between density and refractive index is ex- 
pressed under the form: 



Quantitative analysis. The methodology to obtain an estimate of the 
density, viscosity and thickness of the soft layer is based on the principle 
described here after, for a potential value of the film thickness, the acoustic 
attenuation AA and the acoustic phase shift A0 correspond, each of them, 
to a set of values for the density and the viscosity as obtained from the 
model described above. Considered together, the two sets present a common 
point, thus resulting in a couple values [pc (h), rj c (h)] as shown in Fig. [7| 



Pc ~ Pf 



n c - n F 



(8) 



Pl ~ Pf 



n L - n F 
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where the lines of expected amplitude and phase shifts are plotted at a given 
thickness for a given value of the amplitude and phase shifts versus density 
and viscosity. For a set of values for the thickness, a curve of possible values of 
viscosity and density are obtained. An unique solution [h, pc, rjc] is obtained 
thanks to the SPR data: for each value of thickness at a given SPR angle shift 
A9 corresponds an unique value of the refractive index [3H] that is further 
related to the density of the layer thanks to Eq. |H1 which assumes an identical 
evolution of the refractive index and of the density of the composite layer. 
The methodology is schematically depicted in the following diagram: 



A0 
AA 



[pc(h),r) C (h)} 



A6 

Our method has the unique characteristics to provide the thickness, the den- 
sity and the viscosity of the adsorbed layer out of three independent ex- 
perimental measurements (phase and insertion loss of the SAW and coupling 
angle of the SPR), under assumed values for the density pi and the refractive 
index of the component and under the further assumption of an entirely 
viscous interaction between the composite layer and the acoustic wave (i.e. 
uric > He)- 
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3 Experimental section 

Acoustic waveguide sensor. We used a Love mode SAW sensor obtained 
by coating a 0.5 mm thick ST-cut quartz substrate with a 1.2 fim thick layer 
of silicon dioxide by plasma enhanced chemical vapor deposition (Plasmalab 
100, Oxford Plasma Technology, UK). The interdigital transducers (IDTs) 
were etched in a 200 nm thick aluminum layer sputtered on the substrate 
consisting of 100 pairs of split fingers with a periodicity of 40 /im and an 
acoustic aperture of 3.2 mm. The acoustic filter had its central frequency 
at 123.5 MHz. The sensing area between the transducers was covered with 
a 10 nm titanium, 50 nm gold coating that was electrically grounded to 
shield the acoustic waveguide. The sensor was packaged using a technique 
described elsewhere |37] that defined a 180 to 200 /d open well above the 
sensing area of D = 4.7 mm length and isolated the sensing area from the 
electrode area in order to reduce electrical interactions between the IDTs and 
the liquid. The open well formed a static liquid cell in which the solution 
was manually injected using a micropipette. The SAW phase and inser- 
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tion loss were continuously monitored in an open loop configuration with an 
HP4396A Network Analyzer and logged on a PC through the GPIB interface. 

Combined SAW/SPR setup. The setup has been described elsewhere 
36J. Briefly, a commercial Ibis II SPR instrument (IBIS Technologies BV, 
Hengelo, The Netherlands) based on a 670 nm laser diode has been modified 
in order to replace the usual gold-coated glass slide with the packaged Love 
mode SAW sensor. 

Temperature monitoring. A type K thermocouple (Chromel-Alumel) 
junction was made by wire-bonding two 100 /im-diameter wires (Filotex, 
France) and was monitored by an AD595AD chip connected to an HP3478A 
multimeter. The temperature of the solution right above the sensing area 
was continuously monitored during the experiment and synchronously with 
the SAW and SPR measurements. 

Surface preparation. The device was cleaned in UV-O3 for 30 minutes 
and followed by a 3 hours contact coating with octadecanethiol hydrophobic 
self-assembled monolayer obtained from Sigma-Aldrich. 
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poly(N-isopropyl acrylamide) (PNIPAAm) from Polysciences Inc. (War- 
rington, PA, USA) was obtained with an average molecular weight M w = 
40 kD and a polydispersity index (M w /M n ) of 2 as characterized by the 
manufacturer. A solution of 500 ppm polymer in deionized (DI) water (clean 
room grade, 18 Mfi.cm) was prepared, dialyzed to cut-off molecular weight 
of 12-14 kD (dialysis tubing DTV1200 from Medicell International Ltd, UK) 
over night and stored at 4 °C. 

Experimental method. The temperature dependence of the thiol coated 
SAW sensor was first calibrated in DI water by cycling from 20 °C to 40 °C 
(rising temperature rate: 85 ± 2 K/h, decreasing rate: 49 ± 2 K/h). After 
calibration, the 500 ppm polymer solution was injected while the tempera- 
ture was set at 29 ± 1 °C. Following the procedure described by Plunkett 
et al. |20j . the sensing area was washed with twice the volume of DI water 
after 3 hours at a constant temperature of 29 ± 1 °C. A third volume of DI 
water was then injected for the measurement. The temperature was then 
cycled several times from ~ 20 °C to ~ 40 °C to monitor the influence of the 
temperature induced conformation changes of the polymer on the acoustic 
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and optical signals. 



4 Results 

The experimental results of SAW and SPR data are given after subtraction of 
the temperature calibration in Fig. El that provides the thermal cycle for the 
SAW insertion loss I.L., the SAW phase variation A0 and the SPR coupling 
angle variation A8 for the thermal cycling during the calibration step and 
after the polymer adsorption. The temperature calibrations were fitted by 
polynomial interpolations and subtracted from the raw data with expressions 
as a function of the temperature T, in °C, as following: 

• SAW insertion loss I.L. (dB) = 1.27- 1(T 4 T 3 - 1.29- 1 (T 2 T 2 + 0.422T - 
25.8; 

• SAW phase (f) (°) = 2.57T + 45.1 (corresponding to a temperature 
coefficient of frequency - TCF - of 32 ppm/°C); 

• SPR angle shift A9 (m°) = -22. 7T + 1797. 

The observed average shift with temperature is compatible with the tabu- 
lated variation of bulk water refraction index while the average shifts 
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of the SAW data are following more complicated patterns that involve the 
effect of the transducers and the thermal shift of the acoustic propagation, 
not linearly correlated with temperature as in the SPR case. 
The mechanical and optical constants used for the SAW and SPR models 
are given in Table [T] [JTJ [33 HOI S3 The evolution of the acoustic and op- 
tical data per thermal cycle is given in Fig. One temperature cycle starts 
for the lowest temperature and ends at the next temperature minimum. It 
shows the evolution of the ratio of the SAW attenuation to the SAW phase 
(left axis) and the SPR angle A8 (right axis) for the first 4 cycles. 
The SAW/SPR data measured at 5 different points on each cycle (25.0 °C 
and 35.1 °C during heating up, and 34.9 °C, 30.8 °C and 25.0 °C during 
cooling down) are reported in Table 121 the point of 30.8 °C with decreas- 
ing temperature corresponds to the peak value of AA/A0. The simulated 
values obtained with the methodology presented in the theoretical section 
are reported in Table EJ The evolution of the polymer weight concentration 
and the viscosity of the PNIPAAm- water layer is plotted in Figs. ITT1 and 1121 
respectively, function of the layer thickness. 
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5 Discussion 

The ratio of the SAW attenuation shift to the phase angle shift is a good 
indicator of the variation of viscosity or thickness of the film and of confor- 
mational changes that occur in the film, this analysis can be extended to 
films evolving at constant temperature such as during immunoassay reac- 
tions. Taken separately, the increase of the amplitude AA and the phase A0 
indicate that the film is getting more viscous and thicker, while their ratio 
determines which factor is predominant. As given in Fig. El a proportionality 
exists between this ratio and h^fpc / s/Vc] consequently, if density variations 
are small in comparison to fluctuations of the film thickness and viscosity, an 
increasing ratio denotes a decrease of the viscosity or a thickening of the layer 
while a decreasing ratio indicates the opposite, i.e. an increasing viscosity or 
a thinning layer. 

As the film thickness gets over the penetration depth, the SAW signal ra- 
tio saturates and the acoustic wave sees the deposited film as a semi-infinite 
bulk fluid. In such case, viscosity variations are probed through the individ- 
ual variation of phase or insertion loss of the device. Hence, measuring the 
SAW signal ratio for viscous fluids is meaningless since its value is constant 
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and equals a saturation value specific to each particular SAW sensor. In the 
same order of idea, a direct comparison of the SAW phase shift with the 
one obtained from a standard calibration of the sensor with a bulk viscous 
solution is not adequately fitted to the actual event that takes place on the 
sensing surface during the segregation of a soft layer: at each time that the 
film is smaller than the penetration depth, the film thickness plays a crucial 
role and must be determined independently to estimate correctly the rele- 
vant characteristics of the adsorbed film, as presently for the PNIPAAm case. 

PNIPAAm presents a conformational transition temperature, also called 
Lower Critical Solution Temperature (LCST), that identifies a clear sepa- 
ration of the physical characteristics of the polymer as seen qualitatively 
by acoustic and optical techniques. The transition corresponds to a switch 
from hydrophilic to hydrophobic because of the unshielded isopropyl group 
H3 El The LCST noticed by the abrupt change of SAW and SPR 
signals as seen in Fig. El is located at 30.8 ± 0.3 °C for both rising up and 
cooling down. The transition is accompanied by a modification of the solu- 
bility of the polymer in solution and therefore by a gradient of concentration 
from the surface to the volume. Through this paper, it is assumed that the 
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composite layer is PNIPAAm with water as solvent and a step gradient in 
concentrations defines a clear interface between the adsorbed layer and the 
bulk fluid. Practically, it is more likely that the layer-fluid interface displays 
a smoother gradient due to dangling polymer strands and polymer diffu- 
sion. According to the same principle, surface roughness effects known to 
influence the acoustic signal are disregarded at the moment as the surface 
roughness has been measured with a value in the order of 3-5 nm pp j46j . which 
is much smaller than the acoustic penetration depth and the film thickness. 
These specific effects, gradient of concentration and surface roughness, were 
not considered since the presented methodology is not complete enough to 
address such interfaces yet. To that end, further assumptions on the con- 
centration distribution as a function of the depth are needed and must be 
implemented in the model as, for instance, a multilayered composite layer 
with varying mechanical characteristics following a diffusion law. This type 
of implementation requires additional theoretical and experimental work to 
establish and measure concentration gradients and surface roughness in order 
to determine their exact incidence on the quantitative results presented here. 
Such developments are far beyond the scope of this paper. 



28 



Qualitatively and quantitatively, several variations are detected for each 
thermal cycle, with consistency from cycle to cycle. Variations occur in all 
observed parameters - density, viscosity and film thickness - through the 
cycle. The discussion is facilitated by distinguishing 3 regions per thermal 
cycle: (I) from 25 °C to the LCST; (II) from LCST with temperature in- 
creasing to LCST with temperature decreasing; (III) from LCST to 25 °C. 
Fig. El is a pictorial representation of the assumed mechanisms probed by 
the combined SAW/SPR technique and discussed here after. 
In region (I), a 4 wt % PNIPAAm hydrophilic layer absorbed on the surface 
is at equilibrium with the bulk solution. As the LCST is reached, the con- 
formational transition of the polymer results, in region (II), in a radical and 
direct modification of the concentration that jumps almost immediately to a 
denser layer (12 wt % ), thus a high viscosity resulting from stronger chains 
interactions. During the same time, the polymer in solution has a lower sol- 
ubility and is segregating to the layer that continuously grows till a novel 
equilibrium is reached between surface and bulk concentrations. Meanwhile 
and as seen in Fig. the viscosity increases as more chains interactions 
are taking place. This explanation is correlated with the SPR signal that 
presents a sudden variation after the LCST is reached, reflecting the confor- 



29 



mational change of the layer. On the other hand, the abrupt transition is 
seen in the SAW ratio that presents a slow ramp up in region (I) followed 
by a slight decrease in region (II). As mentionned before, a decreasing ratio 
AA/A<p indicates that the viscosity increases or the thickness decreases, the 
situation that prevails here is an increasing viscosity. 

The situation that occurs during cooling is different as observed at the 
transition between regions (II) and (III). As the temperature reaches again 
the LOST, the polymer is back to a hydrophilic behavior and start to swell 
in water. Because of the swelling, the layer thickness slightly increases but 
the PNIPAAm concentration and the viscosity are dropping fast. The effect 
is seen in the SAW signal ratio that is ramping up since a decrease of viscos- 
ity increases the ratio. Immediately after, in the early stage of region (III), 
the layer is effectively swelling and gets thinner while the viscosity remains 
constant, therefore the ratio is coming back to the initial value it presented 
in region (I). The estimation of the viscosity of organic monolayers is a spe- 
cific feature displayed by the combined technique while few methods are able 
to tackle this type of measurement |17j, so an exact comparison of our re- 
sults with literature values is not directly possible although the estimated 
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values determined by simulations are reasonable and of the same order of 
magnitude as the ones reported by Milewska et al. [IB] fo r bulk solutions 
of PNIPAAm although they characterized PNIPAAm solutions with a larger 
molecular weight of 525 kD. 

The asymmetry seen in the intensity of the SAW signal ratio peaks at 
LCST is possibly linked to the large difference in rising and decreasing tem- 
perature rates. It has been reported by Boutris et al. [49 from optical density 
measurements that the scanning rate influences the temperature transition 
because the system is not at equilibrium for large scanning rates; in our 
measurement, the scanning rate is moderate and the LCST is seen at the 
same temperature for both directions of the thermal cycling, indicating that 
the system is at equilibrium at any moment of the cycle otherwise we would 
observe a shift in the location of the LCST as a function of the tempera- 
ture direction. Boutris et al. equally reported that the dissolution of the 
PNIPAAm during cooling is a much slower process for an identical temper- 
ature scanning rate than the abrupt transition noticed during the heating 
up, which might account for the intensity difference of the SAW signal ratio 
seen upon heating up and cooling down. Finally, as the number of cycles 
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increases, a part of the polymer is diffusing away from the surface and will 
not re-enter the surface-volume equilibrium process of the following thermal 
cycle, explaining for the cycle degradation. 

The temperature change might also cause modification to the interpre- 
tation of the data and are important for the calibration. The SPR sees the 
refractive index modification quite clearly and is linearly dependent upon the 
change with a given slope. The SAW situation differs from the SPR: with 
temperature, the water viscosity drops from ~1 cP to 0.6 cP from 20 °C 
to 35 °C as tabulated jSH], corresponding theoretically to lower value of the 
insertion loss and a phase shift but in the same time the propagation char- 
acteristics of the Love mode device are modified by the temperature shift, 
the net result is the 3 rd order polynomial expression for the temperature on 
the insertion loss of the device while it is included in a linear relation for the 
phase. The effects compensate for each other. 

Our observations are in agreement with the results reported by Plunkett 
et al. of a dense wetting layer that grows in thickness above the LCST. Our 
method provides an additional quantitative evaluation of the layer viscosity 
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below and above the transition temperature, and quantitative information 
determining the amount of solvent and the thickness of the adsorbed soft 
layer. 

6 Conclusions 

The combination of surface acoustic wave and surface plasmon resonance 
measurements delivers useful information to estimate the in situ properties 
of surface adsorbed layers. A general method of data extraction has been 
presented that can be applied to layers adsorbed from liquid environments, 
such as in biosensors. The experimental ratio of amplitude shift to phase 
shift of the SAW signal correlated to simulations is an indicator of a rigid or 
a viscous interaction between the layer and the acoustic wave. A soft layer 
is characterized by a viscous behavior that results in a high damping of the 
acoustic wave, thus to a high ratio of the attenuation to the phase shift. The 
ratio is not a sufficient criterion to discriminate between soft and rigid layers 
since for very thin adsorbed layers, the two possible behaviors have a similar 
effect on the acoustic signal. Without data analysis, the experimental ra- 
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tio is an indicator of the evolution of the thickness or of the viscosity of the 
layer: an increasing ratio is related to a growing film or to a viscosity decrease, 
while a decreasing ratio is related to a thinning film or to a viscosity increase. 

With simulations of an equivalent viscoelastic transmission lines model, 
SAW and SPR data are used to assess quantitatively the density, the viscos- 
ity and the thickness of soft layers. The potential of the method has been 
demonstrated with the example of PNIPAAm, a thermally sensitive poly- 
mer. The record of SAW and SPR data during a thermal cycle clearly shows 
a conformational transition at the LCST point measured at 30.8 °C. The 
LCST transition is seen in the SAW signal ratio by a peak value and in the 
same time by an abrupt variation of the SPR signal. Graphs showing the 
evolution of the PNIPAAm weight concentration and viscosity versus the 
layer thickness were built that indicate first a thickening of the layer followed 
by an increase of viscosity then the collapse of the layer. We also conclude 
from this study that the analysis of the adsorbed layer of PNIPAAm can be 
considered under the assumption of a soft layer as determined by the high 
value of the experimental SAW signal ratio ratio which is in agree- 

ment with a viscous interaction between the layer and the acoustic wave. 
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Simulation results validate the soft layer hypothesis as seen by a high water 
content in the layer (between 88 wt % and 96 wt %). 
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Fig. n Schematic of the PNIPAAm structure for a degree of polymeriza- 
tion n. 

Fig. |21 Schematic cross-section of the combined SAW/SPR setup: a SAW 
device is attached with optical index matching oil to a prism for the coupling 
of a LASER light with the gold coated surface of the acoustic device to gen- 
erate the surface plasmon resonance. The arrows indicate the propagation 
of the acoustic wave (dashed line) and the collimated LASER beam (dot- 
dashed line). 

Fig. El Differential element of the viscoelastic mechanical transmission 
line model: the shear stress (T) and the shear wave particle velocity (v) 
propagate through a symmetric differential element of length dr made of an 
inductance L shunted by the series connection of a capacitance C and a con- 
ductance G. 

Fig. IH Schematic representation of the physical layered structure of the 
Love mode SAW sensor. A guided acoustic shear-acoustic wave propagates 
in the x direction in a stack composed of the semi-infinite piezoelectric sub- 
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strate (S), a guiding layer (G) of thickness t, the composite layer (C) of 
thickness h that segregates from a semi- infinite fluid (F). The sensing area 
of the device is at the interface between the guiding layer and the adsorbed 
composite layer. The plane T is a reference plane used to determine the 
propagation characteristics of the guided acoustic wave. For clarity only 
one guiding layer is represented although the practical device is composed of 
several guiding layers. 

Fig. Equivalent transmission line model of the physical structure in the 
direction of resonance (z) given in Fig. 0] The equivalent transmission lines 
of each part of the device implement the model of Fig. |3] and are schemati- 
cally represented by stub lines with the indication of the wavenumber kj Z and 
the characteristic impedance Zj Z . The lines are connected in series to satisfy 
the boundary conditions and the two semi-infinite layers are loaded with a 
matching impedance. Resonance solutions of this network are solution for a 
guided Love mode. 

Fig. EJ Simulated values for the ratio of the SAW attenuation shift to 
the SAW phase shift AA/A<j) as a function of the composite film thickness 

43 



normalized to the penetration depth in this film 5 = \f2r] C / (pc^)- 

Fig. [7| Plot for a given film thickness of the values of normalized density 
and normalized viscosity for a given and constant value of the SAW phase 
shift A0 (solid line) and a constant value of the SAW insertion loss shift 
A A (dashed line). The values of density are normalized to obtain the direct 
correlation with SPR data. The intersection point of the two curves gives for 
that thickness the density and viscosity of the film. The procedure must be 
repeated for all possible values of the film thickness. 

Fig. |SJ Raw SAW and SPR data temperature calibration and tempera- 
ture cycling effect as a function of the temperature for 5 cycles after removal 
of the temperature effect. The arrows indicate the direction of the temper- 
ature cycling. Top graph: SAW insertion loss, middle graph: SAW phase, 
and bottom graph: SPR angle. 

Fig. El Evolution of the SAW attenuation shift to phase shift ratio (dots, 
left axis) and SPR shift (squares, right axis) for the 4 first temperature cycles. 
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Fig. El Pictorial representation of the mechanisms probed by the SAW/SPR. 
Below LCST, the layer contains a large amount of water and because of the 
solubility of the physisorbed PNIPAAm, there is an equilibrium between sur- 
face and bulk concentrations. As the temperature increases above LCST, the 
conformation changes and results in an hydrophobic material collapsing on 
the surface, in addition to a segregation: as a result, a denser and thicker 
layer builds above the surface. As diffusion drives part of the material away 
from the surface, SAW and SPR signals are getting different as the number 
of thermal cycles increases. 

Fig. ^2 Summarizing graph of the evolution of the estimated PNIPAAm 
weight concentration in the adsorbed layer and the thickness for the 5 probed 
temperatures. The background lines are the measured SPR angle shifts. 

Fig. El Summarizing graph of the evolution of the estimated viscosity of 
the adsorbed layer and the thickness for the 5 probed temperatures. 

Table HJ Mechanical and optical constants used for the SAW and SPR 
models for the materials entering the experimental set-up and for a device 
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operating at a frequency of 123.5 MHz. 

Table El Experimental mean values observed during the 4 first tempera- 
ture cycles probed at 5 different points for the thermocouple, the SAW and 
the SPR. The arrows indicate the direction of the temperature cycle. 

Table El Evaluation of the density p, viscosity rj, thickness h, penetration 
depth in the film 5 simulated with the results indicated in Table El The sim- 
ulated thickness to penetration depth ratio h/5 and the experimental SAW 
amplitude shift to SAW phase shift AA/A(p are given for a comparison with 
Fig. 
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Figure 1: Schematic of the PNIPAAm structure for a degree of polymeriza- 
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Material 


Density 


rii i • nr* 

Shear stillness 


Viscosity 


Ac. Velocity 


Refractive 




p (g/cm J ) 


fi (GraJ 


r) (cP) 


V (m/sj 


index n 


Quartz 


2.648 


66.3 





5004 


1.518 


Si0 2 


2.61 


8.31 





1784 


1.5 


Ti 


4.51 


46.7 





4230 


2.76 + 3.84i 


Au 


19.3 


28.5 





1215 


0.14 + 3.6971 


Water 


1 





1 


19.7 + 19.7i 


1.33 


PNIPAAm 


1.4 








1.47 



Table 1: 
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Thermocouple 

Temp. (°C) 


SA 

AA (dB) 




SPR 

A6 (m°) 


25.0 ±0.1 \ 
25.0 ±0.1 / 


-0.17 ±0.04 
-0.14 ±0.02 


-3.60 ±0.38 
-2.72 ±0.74 


146 ± 17 
142 ± 13 


30.8 ±0.3 \ 


-1.05 ±0.43 


-10.0 ±2.4 


501 ± 21 


34.9 ±0.3 \ 
35.1 ±0.1 / 


-1.35 ±0.22 
-0.66 ± 0.45 


-24.0 ±5.2 
-13.2 ±5.4 


902 ± 52 
681 ± 127 



Table 2: 
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Temp. 


PC 


Vc 


h 


5 


h/8 


AA/Acf) 


( ^) 


(g/cm 3 ) 


(cF) 


(nm) 


(nm) 




(ATI /-r-r,A\ 

{dD/iad ) 


25.0 \ 


1.016 


1.36 


20.5 ± 1.3 


59 


0.35 


2.6 ±0.6 


25.0 / 


1.015 


1.25 


21 ± 1.3 


56 


0.37 


3.0 ±0.5 


30.8 \ 


1.022 


1.83 


59 ±3.6 


68 


0.87 


5.9 ± 1.3 


34.9 \ 


1.045 


3.81 


50 ±5.0 


97 


0.51 


3.2 ±0.3 


35.1 / 


1.048 


2.46 


33 ±2.0 


78 


0.42 


2.8 ±0.9 



Table 3: 
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